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Abstract. - Entanglement of spins is analyzed for two electrons extracted from a mixed many 
electron state by projecting onto the two-electron subspace. The concurrence formulae are ex- 
pressed in a compact form for states with a well defined square of the total spin projection. 
As an example, the thermal entanglement for a qubit pair with an anisotropic Heisenberg and 
the Dzyaloshinskii-Moriya interactions in an inhomogeneous magnetic field is given analytically. 
Remarkably, the concurrence of a pair of electrons with antiparallel spins and in a delocalised 
orbital state is given by the scalar product of the state with its spin-flipped state and not with 
the time-reversed state. 



Quantum entanglement is considered a key resource 
for quantum cryptography and quantum computation [1] . 
Quantifying entanglement [2,3] and identifying maximally 
entangled states is thus important in the planning of de- 
vices. On the other hand, the research of the entangle- 
ment and its decoherence can also lead to a better under- 
standing of the foundations of physics, for example of the 
quantum to classical crossover [4] and of the origins of the 
thermodynamic laws [5]. 

Quantum dot arrays, controlled by electrical gating, are, 
due to their scalability, promising candidates for opera- 
tional devices [6,7]. The building blocks are coupled few 
electron quantum dots enabling full control over individ- 
ual electrons [8,9]. In a quantum dot the qubit is usu- 
ally represented by the spin of an electron. Alternatively, 
charge pseudo-spin entanglement in double quantum dots 
can also be exploited [10], but is prone to the decoherence 
due to the Coulomb interaction with the environment. 

The description of electrons by the spin degrees of free- 
dom only is a simplification valid when the electrons are lo- 
calized and the charge fluctuations are negligible. In gen- 
eral, both orbital and spin degrees of freedom are present, 
but if one is interested in the spin entanglement only, one 
should trace out the spatial dependence. Typical exam- 
ples are a recently proposed route to generation of per- 
fectly entangled electron pairs by the use of acoustic waves 



in the surface of a GaAs/AlGaAs structure [11] or elastic 
scattering of electrons in semiconducting carbon nanotube 
structures with orbital degeneracy [12]. For a special case, 
where there are precisely two electrons in a pure state on 
the lattice, the entanglement can be given by simple for- 
mulae expressed in terms of the wave-function [13]. Such 
formulae can be applied to the determination and opti- 
mization of entanglement generation between static and 
flying qubits in one dimensional systems [14, 15]. They 
also enable the analysis of the entanglement between qubit 
pairs in various double quantum dot configurations cou- 
pled to external leads [16]. 

Here we are interested in the entanglement of an elec- 
tron pair extracted from a many-body state, which can 
be, for example, an open system of interacting electrons 
in a solid state structure of several coupled quantum dots. 
In particular, we take that the measurement apparatus 
extracts precisely one electron from each of two non- 
overlapping regions of the structure - domains A and B. 
The state of the system is arbitrary and includes fluctu- 
ations of electrons between the domains or between the 
domains and the environment, which introduces spin and 
charge fluctuations to the subsystem AU B. 

In this letter we express the reduced density matrix of 
two spin-qubits in terms of projected spin-spin correlators 
which allows the analysis of entanglement of qubit pairs 
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extracted from a general many electron state. The corre- 
sponding concurrence is then given explicitly for systems 
conserving the square of the total spin projection, which 
is illustrated by several examples. 

First consider two separated electrons, one from domain 
A and the other from domain B, with spin states labeled 
by s = ±i and t = ±i, respectively. Let the electrons 
be in a pure state expressed in the standard basis j/z) = 
\st) G {ITT) , ITi) , UT) , Ui>} for e {i, 2, 3, 4}, as 



= ^a^|/i). 



(1) 



Because the electrons are in a state completely determined 
by the spin degrees of freedom only, the entanglement 
can be quantified with the entanglement of formation or, 
equivalently, with the concurrence C = 2\axai — 02^3! 
[17]. In general, two spins may be a subsystem of a larger 
system with many degrees of freedom and the subsystem is 
then described by a reduced 4x4 density matrix p. In this 
case the concurrence is given by the Wootters formula [18], 



C = max(0, 2\max - 



(2) 



where \j are the square roots of the eigenvalues of the non- 
Hermitian matrix pp among which \max is the largest, and 
p is the time-reversed density matrix p. 

In this paper we consider the electrons in domains A 
and B as a subsystem of a total system described by a 
density matrix 



Ptot = ^Pn\n){n\. 



(3) 



Then a projective measurement is performed by an appa- 
ratus which extracts an electron pair: one electron from 
A and another one from B in such a way that after the 
projection, the system is in one of the normalized states 
'P\n) I (n\'P\n), where the projector V projects onto the 
subspace where in each of the domains A and B there is 
exactly one electron. What we are left with is the density 
matrix 

V\n){n\V VptotV 



{n\V\n) TiVpto 



(4) 



where qn = P{n\V) is the conditional probability that af- 
ter the measurement the system will be in the projected 
state 'Pjn), i.e., that the electrons were extracted from 
the state \n). This probability is given by the Bayes' 
formula [1], P(n|'P) = P{V\n)pn/ (V) , where P(V\n) = 
{n\P\n) is the conditional probability for a single occu- 
pancy of each of the domains for a particular state jn) 
and (P) = 5^„Pn("|'P|?^) = TrPptot is the probability 
that the apparatus will click, i.e., that the desired two 
electrons will be extracted. 

The projected states P\n) read 



ij,st 



(5) 



where the operators cj^ and ct^ create electrons with spin 
s at sites i ^ A and with spin t at sites j G B, respectively, 
and, being ordinary electron creation operators, obey the 
fermionic rule Cisc|,,, + c\,^,Cis = Su'Sss'- 

The number of sites within the domains is arbitrary. 
States \^n,ij} normalized and represent empty do- 
mains A and B with the rest of the system in an arbitrary 
configuration. In general, these vacuum states may be dif- 
ferent for each of the states \n) and also for each particular 
occupation of pairs of sites within the domains. The 
projector P removes from \n) all components except those 
where each of the domains is occupied by precisely one 
electron and may be written explicitly as 



P = 



Na 

n 



Nb 

n 



TlB 



(6) 



where 



nA{i 



E 



leA(B),s 



Ci cis is the number operator 



for domains A{B) and Na(b) is the maximum possible 
number of electrons in A{B). 

Being interested in the spin entanglement we consider 
the reduced density matrix where only the spin degrees 
of freedom are retained, p = J^^u Pi^'^\f^){^\' with p/^^, = 
P{st){s't') and 



Pist)( 



1 



l.s't' ■ 



(7) 



This formula is useful if the wave functions are known. 
However, in some cases it is possible to determine vari- 
ous correlation functions of the system without an explicit 
knowledge of the wave functions. Then it is advantageous 
to express the density matrix in terms of spin correla- 
tors [19], 



iP) 



(PlSs) 
(pIPb) 
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(SaPI) 
(SaS^) 

(PiPl) 
{Pis^} 
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(Pa^b) 

{PiPi) I 



where {P) = ^gtiPXPs) ^he probability that in the 
subsystem AuB there will be precisely two electrons, one 
in each of the domains. The correlators are expressed as 
the expectation values of projected operators in the sense 
(O) = J2nPn{''^\T^^'P\''^) where O consists oi A-B pairs of 
operators 

qx,y,z _ 1 x,y,z J 

leA{B),ss' 

ieA{B) 

Here a^'^^'^ are the Pauli matrices, S^^^^ — S^^^^iiS"^^^^ 

and ni^s = cj^cis is the electron number operator. 

The evaluation of the correlators is simplified if each of 
the domains A and B consists of one site only, in which 
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case VOV = O, i.e., the states corresponding to empty 
or doubly occupied sites arc projected away by the oper- 
ator O. A prototype example is a Hubbard dimer, i.e., 
two electrons on two sites and described by the Hubbard 
model, as studied by Zanardi [20]. Note, however, that the 
entanglement measures introduced for fermionic systems 
where multiple occupancy is retained [20-24] are different 
from the entanglement of formation studied here. 

The concurrence for the domains A and B is determined 
from the projected density matrix by the Wootters for- 
mula, eq. (2). In general, Xj can be computed numerically, 
but in some cases, due to symmetry the density matrix 
simplifies and analytic evaluation is possible. Such sym- 
metries were exploited in various coupled spin systems on 
a lattice with translational and parity invariance [25-28]. 

In the present case of interacting electrons in coupled 
quantum dot structures, translational and parity invari- 
ance is an exception. Still, an analogous simplification 
is possible in a special case when the density operator 
commutes with the square of the total spin projection for 
AU B, S'' = S% + S^. For such biaxial systems [29] with 

[P,iSr]=0, (8) 

p is a block matrix: pi2 = pis = P24 = P34 = or, 
equivalently, {S'^^'g^) = and {S\(^b)^b(A)) = 0- 

The eigenvalues of pp (which, again, is a block 
matrix) follow trivially from two decoupled blocks cor- 
responding to subspaces with parallel, (5*^)^ = 1 and 
{p. = 1,1^ = 4}, or antiparallel spins, (S^)^ = and {p = 
2,1/ = 3}. The matrix elements p^i, are interrelated [30], 

^ ^Pf^f^Pi^-/, which leads to Xj = ^p^^p„„ ± 
The concurrence is then determined by 

C = max(0,Cn,q|)/^(PlPi), (9) 

St 

Cn = 2\{S\S-^)\~2^{P\pl){PiPi), 

q, = 2\{S\S+)\-2^ {P\P^^){PiPl), 

which represents a generalisation of the result derived for 
the case of precisely two delocalised electrons in a pure 
state [13]. 

For axially symmetric systems, i.e., conserving the total 
spin projection, [p, 5^] = 0, the formula simplifies because 
{S+S+) = and 

C = max(0,Cu)/(P). (10) 

For SU(2) spin symmetric case, [p, Sa + Ss] =0, the con- 
currence is completely determined by a single spin invari- 
ant and C = 2 max (O, - {Sa ■ Sb) /{V) - \). 

In practice, several specific cases are of interest. Let us 
first consider a special case of the total system being in 
a pure state \m) containing only two electrons, i.e., p„ = 
5nm and j^) = |0). We assume the electrons are in a 
state with the amplitudes -0^,^ = ct2fij, '^m,i3 = ^t^fij 



and = OiiXij, ipt^^.j = a4Xij, where (pij and Xij are 
normalized. Then, if ipij = Xiji the concurrence is given 
by C = 2\aiaA - a2a3\/{V) with {V) = J2i, which 
is the pure spin-subsystem result, renormalized due to the 
projection. If ^ijXij — the concurrence is C = 
2 ||aia4| — \a2Ce3\\/{P)- Additionally, if the state |m) is 
an eigenstate of S^, the concurrence simplifies further to 
C^2\a2a3\/{r). 

An interesting case in point is a pure state with a zero 
spin projection, S^\m) = 0, and (P) = 1. Such states 
are important, for example, in the realization of entangled 
flying qubit pairs, when two initially unentangled electrons 
approach each other and the interaction conserves S^. The 
concurrence is given by 2|(m|S'^ 5^1772)1, but can also be 
expressed as 

C = V(J^)2+4(m|SA X Ss|m)2. (11) 

Here (T) = (mjmflip) is the scalar product of the state 
\m) with its spin-flipped state jmflip) = J^\m) where the 
spin-flip operator = -I- S\Sg reverses the spins 

inAuB. 

For a special case (S^ x Sb) = 0, the expression 
eq. (11) resembles the result for a general pure spin-state, 
eq. (1), with the concurrence given by C = |(m|m)|, where 
\m) = T\m) is the time-reverse of |m) [18]. The time re- 
versal operator is given by exp[— i7r(S'^ + Sg)]JC, where 
/C is the complex conjugation operator [31], which for 
the present case of two electrons with antiparallel spins 
gives T = HFKl. For a special case of a pure spin-state, 
where lAmy oc V'm.ij and = p, the conjugation has 
no effect and the concurrence takes the customary form 
C = 2|Q;2a3| = 2|a2a3|- In general, however, if each of 
the domains consists of at least two sites, the overlap of 
the state |m) with the spin-flipped state [mflip) is different 
from its overlap with the time-reversed state \m) because 
the amplitudes in |m) are complex conjugated while in 
[mflip) they are not. To be speciflc, let two electrons be in 
the state 

I*) = li4Al + ^^U^T + AAl + *4iclT)|0>' (12) 

where the sites 1, 2 and 3,4 represent the domains A and 
B, respectively. The correct expression for the concur- 
rence is C = |(\I/|\E'flip)| — 1 while the scalar product with 
the time-reversed state is (^|^) = 0. Thus, in general 
the entanglement of electron pairs is not related to the 
scalar product of a state with its time-reverse but to the 
spin-flipped state only. This apparent disagreement with 
the pure spin result is no paradox; it should simply be a 
warning and a demonstration that the electron pair with 
a general orbital extend can not be described by a pure 
spin-state. 

Finally, we present an example where the biaxial sym- 
metry [p, {S^)^] = arises naturally in a system in thermal 
equilibrium. Consider a pair of electrons in domains A 
and B representing a weakly coupled double quantum dot 
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structure where charge fluctuations between A and B and 
with the rest of the system arc negligible (i.e., {V) = 1). 
Besides the Coulomb interaction, the spin-orbit interac- 
tion may also be present and a general effective Hamil- 
tonian can contain the anisotropic Hcisenbcrg exchange 
interaction, the Dzyaloshinskii-Moriya term [32] and the 
coupling to an inhomogeneous external magnetic field, 



7 QV QV 



rr T QX QX 

n — JxJA'^B 

+Ba • -|- Bb • Sb + 
+D- (Sa X Sb) + Sa-TSb. 



(13) 



We assume the magnetic field and D are parallel to the 
2;-axis and a symmetric tensor T is of the form 



T 



In the basis {|/u)} the Hamiltonian takes the form 




H = 



/ /ill hi4 

/l22 /l23 

/l^3 /l33 

\ /l44 / 



(14) 



which is the most general form of a Hamiltonian com- 
muting with (5'^)^. The corresponding density ma- 
trix describing thermal equilibrium is given hj p = 
exp(— /3_ff)/Tr[cxp(— /3_ff)], where /3 is the inverse temper- 
ature. The matrix elements of the Hamiltonian eq. (14) 
are related to eq. (13) as /iii(44) = Jz/4 ± {Ba + Bb)/2, 
hu = {Jx-Jy)/4 + it/2, /i22(33_) = -J./4±(Ba-Sb)/2, 
/i23 = {Jx + ■^)/4 + iD/2 with Ba{b) = Ba(b)^ and 
D = £)z. 

The problem decouples into two 2x2 subsystems and 
the resulting two pairs of eigenenergies are given analyti- 
cally, 



where the subscript {^v} denotes {14} and {23} for sub- 
spaces with {S^Y = 1 and {S")'^ = 0, respectively. With 
the corresponding sets of eigenvectors, the density matrix 
elements are easily expressed. 



\h^v\ Smh (3X{fj,v} ^-0{haa + huu)/2 



with 



Z = 2e''^(''"+''*^)/2cosh/3;E{i4} 

+ 26-^^^^^+^^^^'"^ cosh /?X{23} . 

The concurrence is given by eq. (9), C = 2max(0, |p23| 

VPnPil, \Pli \ - VP22P33)- 



Some particular cases of this problem, i.e., an isotropic 
Hcisenbcrg model in a magnetic field and spin Hamilto- 
nians with the Dzyaloshinskii-Moriya interaction [33-35], 
were analyzed before. For = Jy and t = the qubit 
pair is axially symmetric and the concurrence simplifies to 
eq. (10) with 

Cn = ( 1/1231 ""^^"^"> e-^^-/^ - e^'^/A /Z. (15) 

In conclusion, we analyzed the spin entanglement of 
electron pairs extracted from a system of electrons in a 
mixed state by projecting it onto the two-electron sub- 
space. The entanglement is quantified by the concurrence 
obtained from the reduced density matrix which is ex- 
pressed in terms of projected spin-spin correlators for the 
measurement domains and normalized by the probability 
that in each of the two measurement domains there is pre- 
cisely one electron. The formalism is appropriate for the 
analysis of the entanglement of formation for the domains 
of open fermionic systems allowing charge fluctuations of 
the subsystems, as is, e.g., a system of coupled quantum 
dots attached to external leads. 

Simplified expressions are derived for systems with a 
good square of spin projection. The result for the most 
general case of the corresponding two qubit system is given 
analytically, which generalizes particular known cases. As 
an example, the thermal state of a double quantum dot 
with anisotropic Heisenberg and Dzyaloshinskii-Moriya in- 
teractions and in an inhomogeneous magnetic field is con- 
sidered. The concurrence is presented in a simple closed 
form. 

Also considered is an electron pair in a pure, but or- 
bitally delocalized state. An example, relevant to the anal- 
ysis of solid state realizations of flying and static qubits, 
demonstrates that the concurrence is given as the overlap 
of the state with its spin-flipped state - but not complex 
conjugated - therefore not by its time-reverse as is the 
case for pure spin-states. 
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